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Abstract 

In a previous paper, the general approach for treatment of algebraic equations of different order in 
gravity theory was exposed, based on the important distinction between covariant and contravariant 
metric tensor components. 

In the present second part of the paper it has been shown that a multivariable cubic algebraic 
equation can also be parametrized by means of complicated, irrational and non-elliptic functions, 
depending on the elliptic Weierstrass function and its derivative. As a model example, the proposed 
before cubic algebraic equation for reparametrization invariance of the gravitational Lagrangian has 
been investigated. 

This is quite different from the standard algebraic geometry approach, where only the parametriza- 
tion of two-dimensional cubic algebraic equations have been considered. Also, the possible applica- 
tions in modern cosmological theories has been commented. 

1 INTRODUCTION 

Recently there has been an increasing interest in finding solutions of the Einstein's equations in terms 
of elliptic and theta functions. In [1] a solution of the Einstein's equations with a perfect fluid energy- 
momentum tensor was found for the special case of the inhomogeneous Szafron - Szekeres cosmological 
model [2, 3,4,5 ] with a metric 

ds * = dt 2 _ e 2a(t,r,y,z) dr 2 _ e 2/9(t,r,»,») (dy 2 + (L1) 

Solutions expressed through ultraclliptic functions for the case of an relativistic gravitational field of 
a rigidly rotating disk of dust have been found in [6, 7, 8]. Further, in [1] important cosmological 

characteristics for observational cosmology such as the Hubble's constant H(t) — -j=f^y and the deceleration 

parameter q = — ■ R ( t ]- R (*) have been expressed in terms of the Jacobi's theta function and of the Weierstrass 
r (t) _ 

elliptic function respectively. Thus solutions in terms of elliptic functions may turn out to be an important 
ingredient in understanding the s. c. "accelerating cosmological models". 
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Since in [1] a solution has been found only for one component of the Einstein's equations and also 
for one additional, limiting case - the FLRW (Friedman - Lemaitre - Robertson - Walker) metric, it is 
important to know whether there are other cosmological metrics, for which solutions in terms of elliptic 
functions exist. 

In [9, 10], instead of seaching new solutions for special cases, a more general approach has been 
proposed. Concretely, it has been shown that the Einstein's vacuum equations and the gravitational 
Lagrangian can be represented in the form of multivariable cubic algebraic equations. In [9] this was 
performed for a specific choice of the contravariant tensor components in the form g lJ = dX l dX° , where 
X 1 = X l (xi 1 x 2l --,x n ) are some generalized coordinates, depending on the initial coordinates xi, X2, --,x. 
In [10], this was performed for a general contravariant metric tensor g v . In both papers, it has been 
assumed that the contravariant metric components differ from the inverse one (i.e. g^g^ — l\ ^ S k ) and 
consequently, a distinction should be made between covariant and contravariant metric tensor components. 
This is the essence of the affine geometry approach, the essence of which have been clarified in the previous 
paper ([10], Part I). 

Most importantly, the derived cubic equations in all cases clearly suggest that in respect to a chosen 
variable, they can be brought to the parametrizable form of the cubic algebraic equation [11, 12, 13] 

y 2 = 4x 3 - g 2 x - g 3 . (1.2) 

Then, in accordance with the standard algebraic geometry prescriptions, one can set up 

x = p(z) ; y = p\z) , (1.3) 

.92 = 60 ^4 and g 3 = 140 Zfi are com pl ex numbers, called the Eisenstein series and p(z) denotes 
ucr ucr 

the Weierstrass elliptic function p(z) = 77 + X) 

to in the complex plane. 

Previously [9], it has been shown how a selected variable dX l of the cubic algebraic equation 

dx'dx 1 (pv\ l9kr dx k - r ik9lr d 2 x k - r l(l9k)r d 2 x k ) - dx l dx l R a = o (i.4) 

of reparametrization invariance of the gravitational Lagrangian and the ratio | of two coefficient func- 
tions in the performed linear - fractional transformation can be parametrized with the Weierstrass function 
and its derivative. 

The main purpose of this second part of the paper will be to demonstrate how this parametrization can 
be extended to all the variables dX l in the cubic algebraic equation and thus a non-plane, multivariable 
cubic algebraic equation will also be possible to be parametrized, but in terms of more complicated 
functions. This is necessary to be performed in order to find how the contravariant metric components 
g l i = dX % dX^ are parametrized - the importance of this will become evident in the next (third) part of 
this paper. 

It will be demonstrated in section 2 how an embedded system of cubic algebraic equations can be 
obtained, and each equation from this sequence will be for the algebraic subvariety of the solutions - i.e. 
the first equation will be for the algebraic variety of n variables, the second one - for (n — 1) variables,...., 
the last one will be only for the dX 1 variable. In section 3 it has been proved that yet the situation 
is different from the standard case - the obtained functions for the differentials dX 1 , dX 2 , dX 3 are 
complicated irrational functions of the Weierstrass function p(z), its derivative p (z), the (covariant) 
metric tensor and affine connection components, which depend on the global coordinates X 1 , X 2 , X 3 . 
Therefore, it might seem that due to the dependence on the global coordinates one cannot assert that the 
functions dX 1 , dX 2 , dX 3 are " uniformization functions", since they should depend only on the complex 
coordinate z. 

In the next (third) part of this paper it shall be proved that yet the uniformization functions for the 
cubic equation can be found since the generalized coordinates X 1 , X 2 , X 3 can be expressed as solutions 



, where the summation is over the poles 
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of a first - order system of nonlinear differential equations and thus the dependence of the differentials 
on X 1 , X 2 , X 3 is eliminated. 

Throughout the whole paper it has been assumed that the second differentials of the generalized 
coordinates are zero d X z = 0. This may seem to be a serious restriction, but it was necessary to be 
imposed in order to be able to construct the algorithm for finding the solutions of the cubic equation. 
The important point is another - since the main tool of the proposed approach is the linear - fractional 
transformation, will this approach be applicable also for the more complicated cases of 1. the cubic 
algebraic equation of reparametrization invariance, but with d 2 X l ^ and 2. the system of cubic and 
quartic Einstein's algebraic equations. In the subsequent parts of this paper, it will be shown that the 
approach continues to work, although with some modifications and technical difficulties 

2 EMBEDDED SEQUENCE OF ALGEBRAIC EQUA- 
TIONS AND FINDING THE SOLUTIONS OF THE 
CUBIC ALGEBRAIC EQUATION 

The purpose of the present subsection will be to describe the method for finding the solution (i. e . the 
algebraic variety of the differentials dX 1 ) of the cubic algebraic equation (1.4) (in the limit d 2 X k = 0). 
The applied method has been proposed first in [9] but here it will be developed further and applied with 
respect to a sequence of algebraic equations with algebraic varieties, which are embedded into the initial 
one. This means that if at first the algorithm is applied with respect to the three-dimensional cubic 
algebraic equation (1.4) and a solution for dX 3 (depending on the Weierstrass function and its derivative 
is found), then the same algorithm will be applied with respect to the two-dimensional cubic algebraic 
equation with variables dX 1 and dX 2 , and finally to the one-dimensional cubic algebraic equation of the 
variable dX 1 only. 

The basic and very simple idea about parametrization of a cubic algebraic equation with the Weier- 
strass function [11-13] can be presented as follows: Let us define the lattice A = {mwi + no> 2 | m, n e Z; 
u)\,u)2 € C, Im^ > 0} and the mapping / : C/A^> CP 2 , which maps the factorized (along the points of 
the lattice A) part of the points on the complex plane into the two dimensional complex projective space 
CP 2 . This means that each point z on the complex plane is mapped onto the point (a;, y) = (p(z), p (z)), 
where x and y belong to the affine curve (1.2). In other words, the functions x = p(z) and y = p (z) 
are uniformization functions for the cubic curve and it can be proved [13] that the only cubic algebraic 
curve (but with number coefficients!) which is parametrized by the uniformization functions x = p{z) 
and y — p (z) is the above mentioned affine curve. 

In the case of the cubic equation (1.4), the aim will be again to bring the equation to the form (1.2) 
and afterwards to make equal each of the coefficient functions to the (numerical) coefficients in (1.2). 

In order to provide a more clear description of the developed method, let us divide it into several 
steps. 

Step 1. The initial cubic algebraic equation (1.4) is presented as a cubic equation with respect to 
the variable dX 3 only 

A 3 (dX 3 ) 3 + B 3 (dX 3 ) 2 + C 3 (dX 3 ) + G^(dX 2 ,dX 1 ,g lJ ,r'y j ,R lk ) = , (2.1) 

where naturally the coefficient functions A 3 , B 3 , C 3 and depend on the variables dX 1 and dX 2 of 
the algebraic subvariety and on the metric tensor gij , the Christoffel connection and the Ricci tensor 
R ■ 

A 3 = 2pT r 33 g 3r ; B 3 = 6pT r a3 g 3r dX a - R 33 , (2.2) 

C 3 = -2R a3 dX a + 2p(T r ap g 3r + 2T r 3f3 g ar )dX a dX f3 . (2.3) 
The Greek indices a, (3 take values a,(3 — 1,2 while the indice r takes values r — 1, 2, 3. 
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Step 2. A linear- fractional transformation 

dX 3 = 



3 _ a 3 (z)dX + b 3 (z) 



c 3 {z)dX + d 3 (z) 



(2.4) 



is performed with the purpose of setting up to zero the coefficient functions in front of the highest (third) 
degree of dX . This will be achieved if G^(dX 2 , dX 1 , 9lj ,1* , R lk ) = where 



Q = A 3 a\ + G 3 c\ + B 3 a 3 c 3 + 2c 3 d 3 C 3 



(2.5) 



which gives a cubic algebraic equation with respect to the two-dimensional algebraic variety of the vari- 
ables dX 1 and dX 2 : 



P T r j(a g fi)r dX^dX a dX p + K^dX a dX^ + K<pdX a + 2p (^j T r 33 g 3r = 
and K*£p and Ka^ are the corresponding quantities [9] 



K 



S = -R af 3 + 2p^(l + 2^)(2T ap g 3r + T^gpr) 

^3 ^3 



and 



C3 



3p-K 3 g 3r - (1 + 2^)R a3 

C3 C 3 



(2.6) 



(2.7) 



(2.8) 



Note that since the linear fractional transformation (with another coefficient functions) will again be 
applied with respect to another cubic equations, everywhere in (2.4 - 2.7) the coefficient functions a 3 (z), 

b 3 (z), c 3 (z) and d 3 (z) bear the indice "3", to distinguish them from the indices in the other linear- 

— 3 

fractional tranformations. In terms of the new variable n 3 = dX the original cubic equation (1.4) 
acquires the form [9] 

(2.9) 



n 



Pi(n) tti + P 2 (n) m z + P 3 (n) m + P A {n) 



where Pi(n) , P2{n), Pz{n) and P\{n) are complicated functions of the ratios ^ and A 3 ,B 3 , C 3 (but 



not of the ratio which is very important). The variable m denotes the ratio ^ and the variable n is 
related to the variable n through the expresssion 

n = \fh 3 \/C 3 

where 



n + L\ > — 

^3 



(3) 



_ b 3 c 3 c 3 
d 3 d 3 d 3 



L< 3 > = 



(3) 



2 f- 



f + 2 



(2.10) 

(2.11) 
(2.12) 



(3) (3) 

The subscript "3" in L\ ' and I/ 2 

means that the corresponding ratios in the R. H. S. also have the 
same subscript. Setting up the coefficient functions Pi(n) , Pi(n), P 3 (n) equal to the number coefficients 
4, 0, —.92, —.93 respectively, one can now parametrize the resulting equation 



n 2 = Am 3 - g 2 m - g 3 



according to the standard prescription 

dp 



n = p(z) = 



dz 



03 
C3 



= m = p(z) 



(2.13) 



(2.14) 
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Taking this into account, representing the linear-fractional transformation (2.4) as (dividing by C3) 



aidX 3 + & 

dx 3 = cs_ cs. (2 15) 

and combining expressions (2.10) for n and (2.15), one can obtain the final formulae for dX 3 as a solution 
of the cubic algebraic equation 

dX 3 = ™ - 1 C ;, 1 ' . (2. 16) 

d 3 , P (z) r (3) B 3 r (3) V ^ 



In order to be more precise, it should be mentioned that the identification of the functions P\{n) , P 2 (n), 
Ps(n) with the number coefficients gives some additional equations [9], which in principle have to be 
taken into account in the solution for dX 3 . This has been investigated to a certain extent in [9], and will 
be continued to be investigated. Here in this paper the main objective will be to show the dependence 
of the solutions on the Weierstrass function and its derivative. Since only the ratios ^ and % enter these 
additional relations, and not 2 (which is related to the Weierstrass function), they do not affect the 
solution with respect to p(z) and p (z). 

Since B 3 and C3 depend on dX 1 and dX 2 , the solution (2.16) for dX 3 shall be called the embedding 
solution for dX 1 and dX 2 . 

Step 3. Let us now consider the two-dimensional cubic equation (2.6). Following the same approach 
and finding the "reduced" cubic algebraic equation for dX 1 only, it shall be proved that the solution for 
dX 2 is the embedding solution for dX 1 . 

For the purpose, let us again write down eq. (2.6) in the form (2.1), singling out the variable dX 2 : 

A 2 (dX 2 ) 3 + B 2 (dX 2 ) 2 + C 2 (dX 2 ) + G {1 \dX\ 9lJ ,T k l3 ,R lk ) = Q , (2.17) 

where the coefficient functions A 2 ,B 2 ,C2 and are the following: 

A 2 = 2 P T r 22 g 2r ; B 2 = Kg } + 2p[2T r 12 g 2r + r 22 g lr ]dX 1 , (2.18) 

C 2 = 2 P [r r 11 g 2r + 2T r 12 g lr )(dX 1 ) 2 + (K$ +K^)dX 1 + K ( 2 2) , (2.19) 
G 1 = 2pT r 11 g lr (dX 1 ) 3 + K^idX 1 ) 2 + K^dX 1 + 2pp 3 (z)T r 33 g 3r . (2.20) 

Note that the starting equation (2.6) has the same structure of the first terms, if one makes the formal 
substitution —R a /3 — > K^l in the second terms, but eq. (2.6) has two more additional terms K^dX 1 + 
2pp 3 (z)T 3393r- Therefore, one might guess how the coefficient functions will look like just by taking into 
account the above substitution and the contributions from the additional terms. Revealing the general 
structure of the coefficient functions might be particularly useful in higher dimensions, when one would 
have a "chain" of cubic algebraic equations. Concretely for the three-dimensional case, investigated here, 

C 2 in (2.19) can be obtained from C3 in (2.3), observing that there will be an additional contribution 

(2) 

from the term K a ' dX a for a = 2. Also, in writing down the coefficient functions in (2.1) it has been 
accounted that as a result of the previous parametrization 2a = p( z ) _ 

Since eq. (2.17) is of the same kind as eq. (2.1), for which we already wrote down the solution, the 
expression for dX 2 will be of the same kind as in formulae (2.16), but with the corresponding functions 
A 2 ,B 2 ,C 2 instead of A 3 ,B 3 ,C3. Taking into account (2.18 - 2.19), the solution for dX 2 can be written 
as follows: 

dx2 = ^p(z)p(z)VC- 2 + hijdXy + h 2 (dXi) + h 3 
■^p'(z)VC^+h(dX^ 2 + l 2 {dX^ + l 3 
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where hi,h 2 ,h 3 ,h,l 2 ,l 3 are expressions, depending on aj,^,^ (r = 1, 2, 3 ; a, (3 = 1,2), g a/3 , K$ , 

K^i and on the Weierstrass function. They will be presented in Appendix A. 

The representation of the solution for dX 2 in the form (2.21) shows that it is an "embedding" solution 
of dX 1 in the sense that it depends on this function. Correspondigly, the solution (2.16) for dX 3 is an 
"embedding" one for the variables dX 1 and dX 2 . 

Step 4. It remains now to investigate the one-dimensional cubic algebraic equation 

A^dX 1 ) 3 + B 1 (dX 1 ) 2 + C 1 (dX 1 ) + G^(g ij ,T^,R ik ) = , (2.22) 

obtained from the two-dimensional cubic algebraic equation (2.17) after applying the linear-fractional 
transformation 

dX 2 = + (2 . 23) 

c 2 (z)dX +d 2 (z) 

and setting up to zero the coefficient function before the highest (third) degree of (dX 2 ) 3 . Taking into 
account that as a result of the previous parametrization |J — p{z) , the coefficient functions Ai,B\, Ciand 
Di are given in a form, not depending on dX 2 and dX 3 : 

A 1 = 2pT r u g lr , (2.24) 
B\ = F 3 p(z) + K$ = 2p(l + 2—)[2T r 12 gi r + T r ll g 2r ]p(z) + K$ , (2.25) 
& = F lP 2 (z) + F 2 p(z) + K[ 2) = 2p[2Y\ 2 g 2r + T T 22 g lr ]p 2 {z) + 

+ {1 + 2 ^)(k£>+k£ ) ) P (z) + k[ 2) , (2.26) 
c 2 

G° = 2p[Y r 22 g 2r + T r 33 g 3r ]p 3 (z) + K^p 2 {z) . (2.27) 

The solution for dX 1 can again be written in the form (2.16), but with ^, ^, L 2 X \ fciand B\, C\ 
instead of these expressions with the indice " 3" . 

Taking into account formulaes (2.24 - 2.27) for A\,B\ and C\, the final expression for dX 1 can be 
written as 

, -^p(z)p'(z)Jf iP 2 + F 2 p(z) + K[ 2) + f lP 3 + f 2P 2 + hp + h , x 

dX 1 = ^ V , (2.28) 

-^p'(z)^F lP 2 (z) + F 2 p(z) + K? ] +g lP 2 (z) +g 2 p{z) + g 3 

where F\,F 2 ,fi, f 2 , f 3 , fi, gi, g 2 and 53 are functions (also to be given in Appendix A), depending on g a p, 
T r nf) (a,/3= 1,2) and on the ratios ^, ^. 



3 A PROOF THAT THE SOLUTIONS dX 1 ,dX 2 AND dX 2, 
ARE NOT ELLIPTIC FUNCTIONS 

Proposition 1 The expressions (2.21) for dX 2 and (2. 28) for dX 1 do not represent elliptic functions. 

Proof: The proof is straightforward and will be based on assuming the contrary. Let us first assume 
that dX 1 is an elliptic function. Then from standard theory of elliptic functions [11] it follows that 
dX 1 (being an elliptic function by assumption) can be represented as 

dX 1 = Kx{p) + p (z)K 2 (p) , (3.1) 
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where K\(p) and K 2 (p) depend on the Weierstrass function only. For convenience one may denote the 
expressions outside the square root in the nominator and denominator as 

Z 1 (p) = f 1 p 3 (z) + f 2 p 2 (z) + f 3 p(z) + f 4 , (3.2) 

Z 2 (p) = gi p 2 (z) + g 2 p(z)+g 3 . (3.3) 
Then, setting up equal the expressions (3.1) and (2.28) for dX 1 , one can express the function K 2 (p) as 

(1 - K 1 (p))pp'\/F 1 p 2 + F 2 p + K{ 2) +Vk[Z 1 (p) - Z 2 (p)K 1 (p) 



P J F lP 2 + F 2 p + k {2) + Vhz 2 ( P ) 



(3.4) 



The R. H. S. of the above expression depends on the derivative p , while the L.H. S. depends on p only. 
Therefore the obtained contradiction is due to the initial assumption that dX 1 is an elliptic function, 
having the representation (3.1). In quite a similar way, it can be proved that dX 2 and dX 3 are not elliptic 
functions. In fact, this is obvious since they are embedding functions of dX 1 , which is not an elliptic 
function. 



4 DISCUSSION 



The most important result in this paper is related to the possibility to find the parametrization functions 
for a multicomponent cubic algebraic surface, again by consequent application of the linear-fractional 
transformation. The parametrization functions represent complicated irrational expressions of the Weier- 
strass function and its derivative, unlike in the standard two - dimensional case, where they are the 
Weierstrass function itself and its first derivative. The advantage of applying the linear- fractional trans- 
formations (2.4) and (2.23) is that by adjusting their coefficient functions (so that the highest - third 
degree in the transformation equation will vanish), the following sequence of plane cubic algebraic equa- 
tions is fulfilled (the analogue of eq.(65) in [9]): 

A (3 V(3)K 3) + P 2 ] ^m) m h + P^(n {3) )m {3) + P<g = , (4.1) 

A (2 V(2)K 2) + P?\n {2) W (2) + P 3 {2 \n (2) )m (2) + Pg = , (4.2) 

P 1 (1) (n (1) )mf 1) + P 2 (1) (n(i)Ki) + P 3 W (n {1) )m {1) + Pg = , (4.3) 

where m.(3), mp), ™(i) denote the ratios — , — , — in the corresponding linear - fractional transformations 

— 3 — 2 — 1 

and ri(3), np), n>(i) are the "new" variables dX , dX , dX . The sequence of plane cubic algebraic 
equations (4.1 - 4.3) should be understood as follows: the first one (4.1) holds if the second one (4.2) is 
fulfilled; the second one (4.2) holds if the third one (4.3) is fulfilled. Of course, in the case of n variables 
(i. e. n component cubic algebraic equation) the generalization is straightforward. Further, since each 
one of the above plane cubic curves can be transformed to the algebraic equation (i = 1, 2, 3) 

nfi) = Pi ] (n (i) )m 3 (i) + P ( 2 (n (i) )m 2 (i) + P ( 3 l) (n (i) )m (i) + (n (i) ) (4.4) 

and subsequently to its parametrizable form, one obtains the solutions of the initial multicomponent 
cubic algebraic equation. 
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5 APPENDIX A: SOME COEFFICIENT FUNCTIONS IN 
THE FINAL SOLUTIONS FOR dX\ dX 2 , dX 3 IN SEC- 
TION 2 



The functions /ii,/i2,/i3 (depending on the Weierstrass function p(z)) and the functions h,l2,fo (not 
depending on p(z)) in the expression (2.23) for the solution dX 2 of the cubic algebraic equation are 

h x = 2pfe- 2^4 2) p(*)) ( 2r i29i, + r r n 92r) , (Al) 



- 2pL™p(z) (2T r 12 g 2r + T r 22 g lr ) , (A2) 
h b 2 (2) \ (2 ) (2) , , (1) 



^3 =\^- 2-±L?>p(z) j K?> LY>p(z)K£ , (A3) 

h = 2 P ( J - 2^4 2) ) (2r^ 2 . 9lr + r; lfl2r ) , (A4) 



— -2- 

C2 d 2 

(2) 



'3 = 



- 2pL^ (2r^.g 2r + r£ 2 </ lr ) , (A5) 
h 2^4 2) Vf } -4 2) ^ • (A6) 



,c 2 d 2 

The functions Fi, F 2 , F 3 , /1, / 2 , fs, fi, g\, g~2, 53 in the solution for dX 1 are the following 

F = 2p (2^ 252r + T r 22 g lr ) ; F 2 = (l + 2^j (k$ + K$>) , (A7) 

F 3 = 2p (\ + 2^ (2F 12fflr + r r ll92r ) , (A8) 

; h^F 2 -L? , (A9) 
d\ ci 

/ 2 ^ ^ - L«F 3 - 2^4^ , (A10) 
J, s f* - 2^ ft - . (A12) 



ci di 

-4 

ci di 



«, S ('*-2^jf-« . (A13) 
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